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ABSTRACT 

Chandra and XMM-Newton have resolved the 2-8 keV X-ray Background (XRB) into point 
sources. Many of the fainter sources are obscured active galactic nuclei (AGN) with column 
densities in the range of 10 22 — 10 23 cm' 2 , some of which have quasar-like luminosities. 
According to obscuration models, the XRB above 8 keV is dominated by emission from 
Compton-thick AGN, with column densities exceeding 1.5 x 10 24 cm~ 2 . Here, we consider 
whether Compton-thick quasars are detectable by Chandra and XMM-Newton by their direct 
(i.e. not scattered) X-ray emission. Detectability is optimized if the objects individually have 
a high luminosity and high redshift, so that the direct emission has a significant flux in the 
observed band. Using a simple galaxy formation model incorporating accreting black holes, 
in which quasars build most of their mass in a Compton-thick manner before expelling the 
obscuring matter, we predict that moderately deep 100 ks Chandra and XMM-Newton expo- 
sures may contain a handful of detectable Compton-thick quasars. Deep Ms or more Chandra 
images should contain 50-100 distant, optically-faint, Compton-thick sources. In passing we 
show that radiation pressure can be as effective in expelling the obscuring gas as quasars 
winds, and yields a black hole mass proportional to the velocity dispersion of the host bulge 
to the fourth power. 
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1 INTRODUCTION 

Most of the 0.5-2 keV soft X-ray Background (XRB) has been re- 
solved by ROSAT (Hasinger et al 1999; Lehmann et al 2001) and 
recently most of the harder 2-7 keV XRB has been resolved by 
Chandra (Mushotzky et al 2000; Giacconi et al 2001; Barger et 
al 2001; Brandt et al 2001) and XMM (Hasinger et al 2001). The 
source identifications are incomplete at the present time but indi- 
cate that about one third of the sources are associated with normal 
quasars, another third with optically-bright galaxies and the last 
third with optically-faint galaxies (Mushotzky et al 2000; Barger 
et al 2000; Alexander et al 2001). X-ray hardness ratios indicate 
that the fainter sources are generally harder and consistent with 
absorbed sources. This confirms obscuration models for the XRB 
(Setti & Woltjer 1989; Madau, Ghisellini & Fabian 1994; Comastri 
et al 1995; Gilli, Salvati & Hasinger 2001). 

Detailed studies, including X-ray spectroscopy of the brighter 
sources (e.g. Crawford et al 2001) have however shown that the 
column densities of absorbing material lie in the range of 10 22 to 
a few times 10 23 cm -2 . Such columns are insufficient to produce 
the vl v peak in the XRB at about 30 keV (Marshall et al 1982). 
In order to match this peak, sources with columns exceeding 1.5 x 
10 24 cm -2 are required, i.e. sources which are Compton thick. 

Here we address this issue and ask whether such sources 
should have been observed in Chandra and XMM images, and 
whether their apparent absence indicates a problem for the obscura- 



tion model of the XRB. Compton-thick sources are certainly found 
in the local Universe, with two of the nearest active galactic nuclei 
(AGN), the Circinus galaxy and NGC4945 (Matt et al 2000, and 
references therein) being in this class. Good examples of powerful 
Compton-thick AGN are NGC6240 (Vignati et al 1999) which is at 
low redshift and IRAS 09104+4109 which is at redshift z = 0.442 
(Franceschini et al 2000; Iwasawa et al 2001). The main issue is 
whether Compton-thick quasars are common. 

Sources which are mildly Compton-thick (10 24 < A^h < 
10 25 cm -2 ) show direct emission peaking above 10 keV. When 
they are heavily Compton thick (A'h > 10 25 cm -2 , this emission 
peak is down-scattered by Compton recoil and is therefore highly 
degraded. How much is seen in the 2-8 keV band from low redshift 
Compton-thick objects depends almost completely on how much 
emission is scattered into our line of sight. For example, the direct 
line of sight to the well-known obscured AGN NGC 1068 has a 
column density exceeding 10 25 cm -2 (Matt et al 1997), the X-ray 
emission seen is all scattered flux. The z = 0.92 hyperluminous 
IRAS galaxy F15307, known to host a quasar from observed scat- 
tered (polarised) optical broad lines, shows no X-ray flux detectable 
by ASCA and ROSAT (Ogasaka et al 1997; Fabian et al 1996). The 
difference between these sources is probably determined by the ge- 
ometry of the obscuring medium, which may be torus-like in the 
case of NGC 1068 and more spherical for F15307. Since synthesis 
models for the XRB spectrum indicate that most (85 per cent) of 
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the AGN power is absorbed (Fabian & Iwasawa 1999), a spherical 
geometry is more likely to be relevant to the sources dominating 
the hard XRB than a torus one. The scattered flux fraction is likely 
to be small. It is the absorbed direct flux which we are interested in 
here. 

Redshift does help to shift so render mildly Compton- 
thick quasars detectable, leading to an inverse K-correction effect 
(Wilman & Fabian 1999). However the sources still need to be lu- 
minous to be detectable. For the spectrum of the XRB, it is the 
mildly Compton-thick objects which dominate at the vl v XRB 
peak. The lowest redshift members of the Compton-thick class con- 
tribute most to the highest energy part of that peak. 

To assess whether Compton-thick quasars ought to be de- 
tected, we have taken one simple obscured AGN model for the 
XRB (from Wilman, Fabian & Nulsen 2000; WFN), and predicted 
the number of Compton-thick sources expected per square degree. 
We have then passed the predicted spectra through the response 
curves of Chandra and XMM to predict the number expected to 
be detected. Our results indicate that some Compton-thick sources 
should be detectable, particularly in the Chandra 1 Ms images 
(Alexander et al 2001; Rosati et al 2001). Chandra does not have 
the collecting area above 8 keV to detect such hard sources read- 
ily; XMM has too much internal background (relevant because of 
its much larger PSF) at those energies. There may yet be renewed 
claims to resolve the XRB when distant Compton-thick sources are 
finally found in reasonable numbers. We argue that our results are 
not strongly dependent on the model adopted. 



2 THE WFN MODEL 

For the sake of completeness, we provide here an outline of the 
WFN model. At its heart is the semi-analytic galaxy formation 
code developed by Nulsen & Fabian (1995), in which the Cole 
& Kaiser (1988) block model is used to simulate the hierachical 
growth of clustering, and the scheme of Nulsen & Fabian (1995) 
is used to treat the behaviour of the gas within a collapsed halo. In 
short, gas within a radius R = Rcf , where the free-fall time is less 
than the cooling time, rapidly forms stars, producing supernovae 
which can then expel some (for normal galaxies) or all (in the case 
of dwarf galaxies) of the remaining gas from the system. Any gas 
at R > Rcf participates in a cooling flow (CF). Newly formed 
normal galaxies are ellipticals unless all the hot gas is able to cool 
before the present or the next hierarchical collapse, in which case 
they are spirals. Any collapse with at most one infalling normal 
galaxy forms a normal galaxy, with the stars of dwarf galaxies con- 
tributing to the spheroid of the new system. Collapses with more 
than one normal galaxy form a group or a cluster, as mergers be- 
tween normal galaxies are ignored. 

Each of the smallest block model units contains a seed black 
hole of mass 1.6 x 10 6 M Q , and when a block collapses all the 
black holes associated with its merging subblocks are assumed to 
merge into a single black hole at the centre of the new galaxy. 
Such nuclear black holes are then fed by Bondi accretion of hot 
gas from the cooling flow atmospheres which surround them, in 
accordance with the model of Nulsen & Fabian (2000). Whilst thus 
accreting, the black hole is assumed to radiate with an efficiency of 
10 per cent, with a 2-10keV luminosity equal to 3 per cent of the 
bolometric one, with a power law plus reflection spectrum. This in- 
trinsic spectrum is then absorbed by the isothermal distribution of 
cold dusty clouds deposited by the cooling flow, as suggested by 
Fabian (1999). The accretion is terminated by wind-driven gas ex- 



pulsion when the mass of the black hole reaches a critical fraction 
of that in the surrounding spheroid, thereby accounting for the ob- 
served correlation between the mass of the remnant black hole and 
its host spheroid. Thereafter the object shines as an optical quasar 
for 9 x 10 7 yr. For full details of the model, see WFN. 

2.1 Wind or radiation pressure? 

The critical mass of the black hole was determined by invoking a 
wind from the central engine (Fabian 1999; see also Silk & Rees 
1998 who use an energy argument rather than the force one used 
here). We note here that the radiation pressure of the absorbed radi- 
ation alone yields a similar limit (i.e. if the force due to the radiation 
~ iabs/c, where the absorbed power L a b s oc Z/Edd)- 

To pursue this in detail, we adopt the scenario of Fabian 
(1999), where a fraction / of the mass in the core of an isother- 
mal galaxy M(< r) — 2v 2 r/G is in cold, X-ray absorbing 
gas. The gas within radius n has been accreted into a black hole 
of mass Mbh = f2v ri/G, leaving a column density N = 
/« 2 /27rGm p ri beyond. Rearranging these formulae means that 

M BH = f 2 v 4 /G 2 imi p N. 

Much of the radiation from the black hole is absorbed by 
the gas beyond n giving rise to a force on the column of 
/iiabs/47rr 2 c, where fi accounts for how much the matter traps 
and reradiates the energy, /i ~ 0.5 for an isolated cloud but may 
be several for the envisaged optically-thick cloud surrounding the 
nucleus. Then at the limit where the outward force due to radiation 
balances the inward gravitational force on the column, 

L abs =47rGMiVm pC /r 1 . 

(The formula assumes that the column is all at one radius; if it ex- 
tends to r max then a factor of log(r max /r) enters into the inward 
force expression.) 

Let L a b s = /il/Edd = fiAmGMm-pCj '<jt , where <jt is the 
Thomson cross section. Then 

N = /1/2/cTT. 

Substituting this into the above black hole mass formula then gives 
Mbh = -L- 

or 

Mbh = 2.3 x 10 9 -fr ( r~ r) * ■ 

hh \200kms- 1 ) 

This is the mass of the black hole when the gas (which fuels it) is 
blown away. The normalization changes slightly, by losing a factor 
of f/ir, if r\ is redefined as the Bondi accretion radius (as in WFN) 
for gas at the galaxy virial temperature. 

It compares well with the results of Ferrarese & Merritt (2000) 
and Gebhardt et al (2000) for reasonable values of /, fi and fa : 

MBH = 1 - 4xl ° 8 (o4)^^" 1 (^bl^) 4M - 

The above result demonstrates that a wind is not essential for 
expelling the gas surrounding a growing black hole. 

In the present model it appears more as if the galaxy po- 
tential determines the black hole mass, rather than the other way 
round. Note however that the gravitational binding energy of 
a galactic bulge, where the velocity dispersion of the bulge is 
300«30okms _1 , is Sbuigc ~ 2x 10~ 6 u 2 00 M b ui ge c 2 . The en- 
ergy from the central black hole -Bagn ~ 5 x 10~ 4 Mb u i gc c 2 . So 
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only one per cent of that energy can have a major effect on the for- 
mation of that bulge and that effect may have occurred when both 
the black hole and galaxy were young. 



3 THE YIELD OF COMPTON-THICK SOURCES 

For the above model with the free parameters set to the values 
given in section 3 of WFN, we present in Figs. 1^1 results on the 
yield of Compton-thick quasars in deep exposures with Chandra 
and XMM. Firstly, Fig. 1 shows as a function of redshift the total 
areal density of Compton-thick AGN in the model with intrinsic 
2-10 keV luminosities in excess of 10 44 ergs -1 . As discussed by 
WFN, however, such sources are obscured by Ah ~ 15At when 
the obscured phase begins (where At = 1.5 x 10 24 cm -2 , the col- 
umn density above which a source becomes Compton-thick), and 
by Ah ~ 3At at its end. Since there is no scattered flux included 
in the model spectra and the direct emission is almost completely 
suppressed for Ah > IOAt (especially below lOkeV), the vast 
majority of the sources contributing to Fig. 1 will not be visible to 
Chandra or XMM. This is borne out by Figs. 2, 3 and 4, described 
below, which were produced by convolving the object spectra in the 
observed frame with 'response functions' giving the effective area 
of the telescope and detector combination under consideration. 

In Fig. 2 we show the areal density as a function of redshift 
of both the Compton-thick and Compton-thin AGN (the latter also 
includes the unobscured quasars), which would produce more than 
30 counts in the 5-10 keV band in a 220 ks exposure with XMM 
on the EPIC pn chip. This is repeated in Fig. 3 for a 100 ks ex- 
posure. Fig. 4 shows the analogous plot for the 2-8 keV band with 
the Chandra ACIS S3 chip, with the count threshold reduced to 
10 counts owing to the smaller PSF (and thus lower internal back- 
ground) compared to XMM. The majority of the detected objects 
contributing to the yields are quite close to the count thresholds. 
For example, if the XMM threshold is increased to 40 counts, the 
yield of Compton-thick (Compton-thin plus unobscured) AGN falls 
over all redshift by a factor 2-3 (1.5-2). Similarly, increasing the 
Chandra detection threshold to 15 counts, results in a fall by a fac- 
tor of 2.5 in the yield of Compton-thick objects, whilst the number 
of Compton-thin and unobscured objects are relatively unchanged. 

Most of the predicted Compton-thick sources should become 
detectable in exposures of lMs (Fig. 5), such as are now being done 
with Chandra (e.g. Alexander et al 2001; Rosati et al 2001). In an 
8.4 x 8.4 arcmin 2 region and 1 Ms exposure we predict 182 sources 
(Compton thick and thin) in the 2-8 keV band whereas Alexander 
et al (2001) find 102. We do not consider that this discrepancy in 
numbers is serious, since most of the sources would be close to the 
detection threshold. That many of the sources are seen in both soft 
and hard energy bands (Alexander et al 2001) is also not a major 
problem since a few per cent of the primary soft radiation may be 
scattered into our line of sight. 

The important issue of the redshift distribution of the sources 
is considered in the next Section, where we also discuss the impli- 
cations of these results for the interpretation of deep Chandra and 
XMM surveys. 



4 DISCUSSION 

We see that the surface density of Compton-thick sources expected 
in a 220 ks exposure (equal in length to the initial Chandra obser- 
vation of the HDF-N and its vicinity by Hornschemeier et al. 2001) 
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Figure 1. The areal density as a function of redshift of all Compton 
thick AGN in the model with intrinsic 2-10 keV luminosity in excess of 
10 44 erg s -1 , i.e. those of quasar luminosity. 



5—10 keV band. 220ksec exposure, sources with > 30 counts 
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Figure 2. The solid line shows the areal density as a function of redshift 
of the Compton-thin and unobscured AGN combined (with no restriction 
on luminosity), showing more than 30 counts in the 5-10 keVband in the 
XMM EPIC pn chip in a 220 ks exposure. The dashed line represents the 
yield of Compton-thick objects. 
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Figure 3. As in Fig. 2 but for an exposure time of 100 ks. 
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2-8 keV band. 220ksec exposure, sources with > 10 counts 
Chandra ACIS S3 
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Figure 4. As in Fig. 2 but for the 2-8 keV band on Chandra ACIS S3 chip, 
and a lower count threshold of 10 counts. 



2-8 keV band. lMsec exposure, sources with > 10 counts 
Chandra ACIS S3 
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Figure 5. As in Fig. 2 but for the 2-8 keV band, a 1 Ms exposure with 
the Chandra ACIS S3 chip, and a lower count threshold of 10 counts. The 
dashed line again shows the yield of Compton-thick objects. 



with Chandra ACIS chip is about 200. A chip is 8 x 8 arcmin 
or 0.0178 sq deg which means that 3-4 faint hard sources are ex- 
pected in such an exposure. If the ACIS-I array is use, with 4 times 
the area, then the number rises to about 10 (given the degradation 
of sensitivity off axis: most of the photons detected are above 6 
keV). The sources have redshifts of 3-7 and are basically detectable 
because of the inverse K-correction. Most of the Compton-thick 
above z — 3 become detectable in deeper exposures of 1 Ms. 

In an XMM pn chip exposure of similar duration, concentrat- 
ing on the 5-10 keV band, we have about 40 sources detected per 
square degree, mostly in the redshift range of 2-5. The field of view 
is about 0.2 sq deg (30 arcmin diameter) leading to about 8 sources 
detected per 220 ks exposure, similar to the yield from Chandra. 
Source confusion should not be a problem at such high photon en- 
ergies. 

The adopted WFN model thus predicts that the direct absorbed 
emission from some Compton-thick quasars is detectable by Chan- 
dra and XMM with similar yields. The absorbed primary emission 
from these sources will be very hard and most will be detected close 
to the threshold for detection (10 counts for Chandra, 30 for XMM), 
as discussed in section 2.1. Scattered emission may however allow 



the sources to be detected also in softer energy bands. The appear- 
ance of such objects in the optical and near infrared bands depends 
on the rate of star formation. If this rate is low then they should ap- 
pear as redshifted early-type bulges, without strong emission lines, 
and if high they will be brighter and of much later spectral type, 
probably with emission lines. 

So far only one source with a spectrum consistent with being 
Compton thick has been reported (Norman et al 2001); its redshift 
z = 3.8. The observed X-ray emission in this case is not the ab- 
sorbed direct radiation but a reflected component. Only a guess can 
be made as to the power of the direct component. It is this direct 
absorbed component of such sources which is important for the 
spectral peak of the XRB. 

Our model predicts that more than half the sources detected in 
a Ms exposure with Chandra will have redshifts ranging from 2-8. 
Most of them would have I band magnitudes exceeding 24, i.e. are 
optically faint according to the definition of Alexander et al (2001). 
The precise magnitudes would depend upon the star formation rate 
in the host galaxy. If we assume that the obscured black hole growth 
phase coincides with the major star formation phase in the galaxy, 
then the star formation rate must be tens Mq yr -1 . Then the R 
(and approximately /) magnitude of the hosts will be about 24 for 
objects at z ~ 3 (e.g. Pettini et al 2001). Objects at z ~ 6 will 
be about 1.5 mag fainter. Alexander et al (2001) have 47 objects 
fainter than I = 24 (15 with / > 25.3), which is about half the 
number we predict. 

At this stage we consider that our model remains viable, with 
a factor of two in number densities arrangeable by a small alter- 
ation of some of the parameters (e.g. the gas fraction and/or the 
gas metallicity: note that most of the sources are predicted to be 
close to the detection threshold, even in lMs). Further work on the 
optically-faint population in the Ms Chandra fields will resolve this 
issue by comparison with the redshift distribution in Fig. 5. 

We now ask how robust our conclusions might be, since they 
depend on a specific model. First we note that they should be a 
reasonable (factor of two) approximation to any model which has 
the obscured phase before the 'normal', unobscured, quasar phase, 
since the model has that emission peaking between z — 1 — 2 
as observed. It also produces a mean local black hole density in 
agreement with that inferred from observation (Merritt & Ferrarese 
2001). The higher redshift (z > 3) Compton-thick sources con- 
tribute 10-20 per cent of the intensity of the XRB spectrum in the 
5-12 keV band, and so are not of great importance to the spectrum; 
the lower redshift sources with 1 < z < 2 contribute far more and 
dominate the spectrum in the 20-30 keV band. These last objects 
are completely undetectable with Chandra or XMM, unless they 
have scattered soft emission. The lowest redshift Compton-thick 
sources contribute most at the highest energies. 

We conclude that distant Compton-thick quasars, which may 
dominate the XRB intensity over the 5-20 keV band and repre- 
sent the growth phase of massive black holes, are just detectable 
with present instruments. A handful may be found in 100 ks expo- 
sures at flux levels close to the detection threshold. 50-100 distant 
Compton-thick sources are predicted to be detectable in Chandra 
deep fields of 1 Ms or more. The sources should be optically faint 
and probably appear as high redshift starburst galaxies. 
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